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Abstract 



(N 
> 

, We made an analysis of the equations of motion which are obtained from the 

one-loop effective action for Brans-Dicke gravity with dilaton-coupled mass- 
less fermions in a time-dependent conformally flat background Various par- 
, ticular solutions, including the well-known stationary one, of the corresponding 

O I set of first-order differential equations are given. Some of these solutions describe 

' an expanding time-dependent Universe with increasing, constant or also decreas- 

ing dilaton. This is illustrated by a numerical analysis. For the nonstationary 
Q^' solutions a stability analysis is given. 
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^ : 1 Introduction 



In recent times Brans-Dicke(BD) gravity with matter in the Einstein frame attracted 
much attention. Especially, such theories have been realized as Einstein gravity with 
dilaton coupled to matter (for a review of the renewed interest in scalar-tensor gravity 
see Ref. 0). Obviously, Brans-Dicke theory |0 represents one of the simplest examples 
of scalar-tensor (or dilatonic) gravities where the background is described by the metric 
and the dilaton. The reason for the consideration of such models is the following: 
First of all, the dilaton is an essential element of string theories, and the low-energy 
string effective action may be considered as some kind of BD theory with higher order 
terms (for a recent review, see [Q). Second, dimensional reduction of Kaluza-Klein 
theories naturally may lead to BD gravity. Third, dilatonic gravity is expected to 
have such important cosmological applications as, e.g., in the case of (hyper) extended 
inflation . In addition, there was some activity on the study of BD cosmologies with 
varying speed of light H . 
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Recently, the anomaly-induced action for dilaton coupled scalars, vectors and spinors 
in four dimensional curved spacetime has been calculated in Ref. 0. Thereafter, in 

the effective action formalism has been applied for a study of quantum cosmology 
in BD gravity with dilaton coupled spinors. There, the one-loop anomaly-induced ef- 
fective action (EA) due to massless fermions on a time-dependent conformally flat 
background coupled to the dilaton has been computed and the (fourth-order) quantum- 
corrected equations of motion have been derived. However, because of the complicated 
structure of that coupled set of (ordinary) differential equations only one special so- 
lution, representing an inflationary Universe with slowly expanding BD dilaton, could 
be presented. 

Here, we study the equations of motion, for conformal as well as physical time, 
more extensively. We were able to show additional particular, physically relevant solu- 
tions. Making use of the symmetry group analysis we showed that the former solution 
is a special case of the "stationary solution" of the set of five first-order differential 
equations being equivalent to the equations of motion. By a numerical analysis the 
asymptotic behaviour of various solutions will be shown; some of them approximate 
the stationary solution. For the non-stationary solutions a perturbative analysis has 
been done and some new solutions corresponding to terminating series are presented. 
Finally, a stability analysis of the equations of motion is given. 



2 Equations of motion: conformal time 



The action (in the Einstein frame) of BD theory with dilaton coupled to A^ massless 
Dirac spinors is 



S 



d^x- 



R 



167rG 



1 ^ 
-{V,<P)iVy) + exp (A0) ^iYV,^' 



i=l 



(2.1) 



with A = — 8A/vrG'/(2co' + 3), uj > —3/2; G is Newtons constant and oo is the BD 
coupling parameter of the dilaton. In the following we restrict ourselves to FRW type 
cosmologies, 

ds'' = -dt^ + a\t)dl\ (2.2) 

where dl'^ is the line metric element of a 3-dimensional space S with constant curvature. 
It is convenient to introduce conformal time rj by means of 



dt = a{ri)dri , 



(2.3) 



to get a space-time which is conformally related to an ultrastatic space-time with 
spatial section M of constant curvature. For the fiat case (S = R^, i.e., k = 0) the 
metric gets 



with 



a{r]) 
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In that case the complete anomaly-induced EA for the dilaton coupled spinor field 
becomes fGl 



W = Vsj dri {haa"" - (61 + b2){a" - (a )')'} 



(2.4) 



where V3 is the (infinite) volume of fiat 3-space, a = o'+A(j)/3 {' = d/dr]) and, for Dirac 
spinors, h = N/10{A'k)'^, 62 = -11A^/I80(47r)2. The total one-loop EA is obtained by 
adding to W the classical action (for A; = 0): 



S=-V3 

2 ^ 



dr] 



A-rG 



'a" + a'2) + 



2o- 



(2.5) 



Then, the corresponding equations of motion are (see Eqs. (15) and (16) of [jT|) 

C^A<t>l?> ^ (3/4^G)a" + a0'2 = , (2.6) 



(A/3)Cae^<^/=^ - (0^0')' = , (2.7) 

where 

a 

with the notations 

~a{j]) = e^^") = a{r]) e^^'') , v{7]) = A<f){7])/3. 

A complete integration of these equations appears to be hopeless. Therefore, let us ask 
for particular solutions. 

First of all, Eq. ( |2.7|) because of ( |2.8|) may be integrated immediately. Then the 
system (p^ ) - (p.8|) , written in terms of the variables d and v, reads: 

(1 - a) [(f " - v'^)d + 2a v'] + aa' = , (2.9) 



S^fi'" - 3Sa a" - 2i/fi'3 + pd^'e-^" = Bh^ , (2.10) 

where B is an arbitrary integration constant, and the following notations have been 
used: 

_ _8 1 ^ ^^_^ = _^>n =^ = _11 

""127rG 3a; + 3/2- ' ^"261^2 aNG - ' ~ h 18' 
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Obviously, a particular solution exists for a = 1 or, equivalently, = 7/6 which, 
because of a" = 0, i.e., a = Ci + C2?7, is given by (a third integration constant for later 



convenience is called c|) 



a(r/) 



CgCi + Csr/ + —{ci + C2r]) + —- ■ 

PC2 13 {ci + C2r]Y 



1/2 



C3 + 



25 



+ 



f3c2 (Ci + 02?]) f3 (Ci + 02?]) 



f2.111 



(2.12) 



This solution describes an expanding Universe with time-dependent dilaton which de- 
creases to some constant value — (3/y4)lnc3; see also Sect. 4 (a2) and Fig. 5. Let us 
now consider the scalar curvature R = 6a"/ of corresponding to that solution. Using 
the following abbreviations 



X 



.1/2, 



ci + C2r]), 



B/{(3c2cf), n = u4/f], 



we obtain 



R 



Yin {n + 4bx^ + 3%^) - b'^/2n 



(2.13) 



C3 (n + 2bx^ + 

In the special case B = 0, which results in 6 = 0, this expression simplifies. However, 
in both cases the asymptotic values for rj ^ oo coincide, 

88 NG 1 



R. 



IStt 3cj + 2 44ri^ ' 



(2.14) 



therefore, depending on the sign of cs, the space of negative (or positive) curvature 
approaches the flat one quite fast. 

Second, since v is not explicitly involved, let us change the variables according to 

s{ri) = \na{ri), wirj) = v'irj)] 

then, instead of Eqs. (p.9|) and ( p.lO| ), we obtain 

s'" + w" - 2{l + u){s' + wf + [3we'^' = B, 

/ I o ' I 2 I I II X i2\ 







(2.15) 



or, equivalently, presupposing a 7^ 0, 1, 

y," = "^(b- f3we^^ + 2uiw + s'f'] 

1 — a V / 

+— f (1 + a)w^ + w' +{2 + ?,a)ws' + 3as 
a V 



a 



w' + 2ws'). 



(2.16) 
(2.17) 



These equations may be subjected to a symmetry group analysis [Q. As a result 
it follows that in the general case {B ^ 0) only the generator of translations along r] 
exists: 

However, in the case i? = an additional generator exists: 

X2 = - wdyj - ds . 

The corresponding invariants are rjw = W and rje^ = g. 

Choosing W = hi = const., g = —h2 = const, we get the particular solution which 
is already known from [|I|: 

w = A^'{ri)/3 = hi/r], (2.18) 

e^('^) = a(r/) = -/i2/r/. (2.19) 

The values of hi and /i2 are determined by 

hl-3hi + 2a = , 
Phj- a) + 2u{hi-l)yhi) = 0, 

leading to 

1 ^ — _^ ^2 1 A 14 ,11 11/ii 

2a! 6a 



(2.20) 



/ii = - (3 ± V9^8^ , /i2 = _(^i-_« + 



from this it is obvious that in order for hi to be real it is necessary that a < 9/8 or, 
equivalently, 47/54 < cu. This solution, corresponding to inflationary universe with 
linearly growing dilaton, has constant scalar curvature R = 12/h2 > 0. 

3 Equations of motion: physical time 

In order to be able to give a physical interpretation of the (particular) solutions let us 
shift back to physical time t: 

dt = a{r]) drj, dr, = a{t)dt . (3.1) 



^Note that the correct expression for Hi = A/3hi, Eq. (22) in Ref. is given by Hi 
-(3/16)V4^G(2^ + 3) (1 T ^^-^1^ 



Then the above Eqs. ( |2.16| ) and ( 2.17 ) change into 



vtt + SstVt + + a{su + 2sf) = , 

sttt + vttt + ^st{stt + Vtt) + (stt + 2sf){st + vt) (3.2) 
- 2{1 + iy){st + vtf + (3vt = Be-^\ 

where now 

s(t) = Ina(t), w{t) = a(t)vt(t) . 
This may be rewritten as follows {a ^ 1): 

^ = St^Vt-U, 

= gft + 2s\ - ^ \ J y^t^ust^ u^) , (^3 3^ 

Bq'^' = utt + "iutst + -^—u{ut + ust + u^) - 2(1 + v)v? + /3(m - St) . 

1 — a 

Obviously, there are 5 constants required to fix any solution at t = to! we chose them 
as follows: 

s(to) = So = Ina(to), St(to) = Si = at(to)/a(to), f (to) = ^^o = ^0(^o)/3, 

M(to) = Mo = Mt(to) + St(to) = t^l + Si, Mt(to) = Ml = Vttit^) + Stt(to) = ^2 + S2. 

(3.4) 

Let us change, by choosing Ut = q and St = p, the above set of differential equations 
into another one being only of first order. Then the basic equations to be studied in 
the following are given by 

Ut = q, St=p, Vt = u-p, 

Pt = -2p^ + - -iq + up + u^), (3.5) 

il-a) 

qt = -Spq-- -{q + up + u'^) + 2{l + u)u^ + p{p - u) + Be-^' . 

(1-a) 

The two constants, go = q(to) and po = p(to), which in addition to Sq, Uq and Vq define 
any solution of Eq. (|3.5| ) are given as ui = go and si = po- 

The variable v is not directly involved into these equations. Hence, one only has 
to find solutions of the four equations for the variables s,u,p,q, and afterwards the 
function v is simply obtained by quadratures. Furthermore, if i? = the same is true 
about s. Therefore, we should distinguish the cases B = and B ^ 0. 
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3.1 Stationary solution 



The set (|3.5|) of first order differential equations for i? = has a stationary solution 
which is of physical relevance. It is determined by 

u = Ust, p = Pst, g = 0, (3.6) 

where Ust and Pst are expressed as follows: 



P» = f- = (3.7) 

Substituting these values of pst, Ust and g = into the equations for s and v, and 
integrating the latter, we obtain the following linear dependencies 

S = So + {t - to)pst , V = Vo+ {t-to){Ust-Pst) ■ (3.8) 

Obviously, this generalizes the solution given by Eqs. (p.l8|) and (p.l9|) , which are 
obtained for to = 0, Sq and Vq arbitrary, pst = 1/^2? ^st = (1 ~ ^i)/^2- The solution 
( |3.8| ) corresponds to an exponentially increasing (pst > 0) universe with a linearly 



growing BD dilaton and constant scalar curvature R = 12pf 



sf 



3.2 Numerical analysis 

In order to get some insight into the various types of behaviour of the solutions of 
the system ( ^.5| ) we made a systematic numerical analysis. Here, we present only the 
general result. 

In the case po > uq, i.e. vi < 0, and B > the solutions for t ^ oo very fast 
approximate the stationary solution (|3.6|) - (|3.8|) . In these cases different values of 
Po,Uo,qo and B change the shape of the solutions only quantitatively; different values 
of N show only a qualitative change for short times but have the same asymptotics; cf. 
Figs. |I] and |^ as well as left panel of Fig. |[ However, if po < uq, i.e., vi > 0, and B < 
the behaviour is quite different, cf. Fig. ^ as well as right panel of Fig. ^. In that cases 
the solution shows eventually (damped) oscillations around the asymptotic solutions 
or exponential increasing (i.e., explosion-type) behaviour; furthermore. Fig. |^ shows a 
dilaton-driven collapse. 

Now, we present some of the plots illustrating the typical behaviour of solutions of 
the system (|3.5|) for different values of parameters and initial values (as in a; = 500 
has been chosen). The dashed line in these plots indicate the stationary solution. 

Having so different types of numerical solutions it is useful to discuss a possibility 
to construct any type of analytical solutions. This will be presented in the next section. 
Also the condition of structural stability is of great importance; it will be considered 
in detail in the last section. 
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Figure 1: Dependence of functions u, p, q, s, v and a = exp(s) for cu = 500, N = 10, 
B — and uq — —1, po — 1, qo — —10. 
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Figure 2: Dependence of functions u, p, q, s, v and a = exp(s) for cu = 500, N = 100, 
B — and uq — 0.1, po — 0.1, = 0.1. 
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Figure 3: Dependence of functions u, p, q, s, v and a = exp(s) for cu = 500, N = 10, 
B — and Uq — 3, po — 2.8, go — 1- 
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Figure 4: Dependence of functions s and v ior u — 500, = 10, Uq = 1, po = ?o = 1 
and different values of B: 

left panel: 1-5 = 0, 2-5 = 10^, 3 - 5 = 10^ 4 - 5 = lO'^. 
right panel: 1-5 = 0, 2-5 = -10, 3 - 5 = -40, 4 - 5 = -60. 
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4 Non-stationary solutions 



In the input equations one can exchange the independent variable t and any of the 
dependent variables. For example, choosing u as new independent variable we can 
rewrite equations ( |3.5| ) in the following form 

qtu = 1 , qSu=P, qvu = u-p, 

qpu = -2p^ + (q + up + u^), (4.1) 

1 — a 

qqu = -^pq-- {q + up + u^) + 2{l + v)u^ + (3 {p - u) + Bq-^'. 

1 — a 

Instead of choosing u as new independent variable we can also choose any other 
dependent variable, for example q. This case is of particular interest because, for g = 0, 
there exists the stationary solution (|3.6| ) of the basic equations ([Ol) . Therefore when 
g 7^ we will get the time-dependent solution according with the first of the following 
equations 

Rtq = 1 , RSq = p , RVq = U — p , 

Rpg = -2p'^ + Y^{q + up + u"^) , Ruq = q, (^4 2) 

u 

where R = -3pq (g + up + u^) + 2(1 + v)u^ + I3{j> - u) + Be'"^^ . 

1 — a 

In view of the polynomial dependencies of ([4. 1|) and (|4.2|) upon u and q it is possible 
to look for solutions in form of an infinite series in these variables. 



4.1 Series in u 

(a) Let us start with the particular case 5 = 0. Then we use the following represen- 
tation for q and p 

00 00 
g = J]A„M", p=Y,BnU^, (4.3) 

n=0 n=0 

while the remaining dependencies are obtained by integrating the first three equations 
in ( fl.lD for t, s and v. Substitution of ( [4.3| ) into ( [4.1[ ) yields an infinite set of equations 
for the coefficients An and Bn'- 

k 

^ (^(/c + 1 - j)Ak+i^jAj + 3Bk-jAj^ - 2(1 + z/)4,3 + PSk,i - (3Bk 

j=0 

+ (Sk,3 + Ak^i{l - 5k,o) + Bk-2{1 - 4,o)(l - = , (4.4) 
1 — a V / 

k ^ 

+ j)Bk+i-jAj + 2Bk-jBj^ - (^4,2 + Ak + Bk-i{l - 4,o)) = . 

j=Q " 
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Of particular interest are those values of the parameters and initial values when these 
series are truncated, i.e., reduce to finite sums. In this case the requirement of vanishing 
for the coefficients Ak and Bk for k > k max yields some additional conditions imposed 
on the coefficients Ak and Bk, that can be fulfilled only for some particular values of 
the parameters involved. Below we present two examples of such solutions, 
(al) The first example is valid for a = and fcmax = 2 and is given by the formulas 

Ao = Bo = Ai = B2 = , A k = Bk = 0, k > 3, 
Ei = l, ^2 = -l±ym^< 0, 

s = So — (IM2) In (1 — A2UQ{t — to)) 5 V = Vq = const. , (4-5) 



q = A2U , u = p 



A2UQ{t - to) 



This solution confirms that different types of behavior of solutions are possible depend- 
ing upon the initial conditions. For mo = po > we have transition to the stationary 
state with zero asjTuptotic values of u and p at t — > 00, while for uq < we have an 
explosion-type behavior whose singularity lies at t — to = l/(^2^^o)- The scalar curva- 
ture is positive and, in the limit t ^ 00, approaches 6 (1 ± \/l + z/)(l =f \/1 -(- 1/)~^ t~^. 

The solution described by formulas ( [4.5|) is valid for a = and corresponds to 
p = u. However, in this case a more general form of the behavior of functions q and t 
upon u (as compared to ( |4.5|) ) can be obtained that is given as follows (C = const.): 

r 

1 



t-to = T^ / dr(v^rTi;-i-r) ' *^(v/rTi;+i+r^ ' 



ro 



So + ^ ^— — In — , v = Vq = const. 

4 VI + y \^f\Tv-\-r; 



q = u^r, u = p = ±C ^Vl + — 1 — r j * "^^'^ ^vT+^ + 1 + r 

(4.6) 

The result of the previous case ( [4.5| ) arises if we impose the restriction = that 
leads to constant values of r = — 1 ± yjl + v. 

Obviously, as long as p = u these solutions are quite special since they correspond 
to a constant dilaton. Furthermore, because oip = u the solution does not depend on 
13 (and G)\ 
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(a2) The second example is valid for a = 1, k max = 4, and is given by the formulas^ 

Ao = Bo = A^ = B2 = A^ = B, = 0, = = , k>5, 
Bi = l, A2 = -2, A, = -B, = 2u/(3, 

111/ u Uo \ 

t — tn = 1 , arctan — , — arctan — , , 

2u 2uo 2,/^\ ^/^j (4.7) 

1, fu\u'-(3/v)\ 1 fu^-l3/v 

" ^ ■^o - 7 In ( 7^ ) , v = Vq + - hi' 



4 V«oK-/5H/ 2 \ul-P/v 

This solution may be obtained also directly from Eqs. ( |2.11| ) and ( |2.12| ) where we 
used the conformal time t]. Indeed, putting i? = in ( p.ll|) and substituting into ( p.3|) 
we get 

dt=^^7W/?) + (ci/ci)g, q = {cr + c,r^)\ 
which is easily integrated as 



t = const. 



■ + I U^/P+{cycl)q - arctan y^-l - (ci/ci)(/?/z/)g 



Making use of the definition oi u = Vt + St = arf/a? + Vr^/a we get, after differentiation 
of ( p.ll| ) and ( |2.12| ), the following relation between u and a: 



u ^ = (ci + C2r])a{r])/c2 = + (ci/ci)g. 

Eliminating q with the help of the last formula from the expression for t we get 



t = const. H — -J —V 1 3 arctan — , 

which coincides with ( |4.7] ) if we take into account arctan x = 7r/2 — arctan(l/,T). The 
rest of Eqs. ( |4.7] ) may be checked immediately. 

In the limit (t — const.) the last formula gives u ^ {t — const.)"^^^ ; in the 

opposite limit t/ \/ —v j P >> 1 we have u ~ (1/2) — const. + {n / A) y^^iTJjfj 

Even in the case of B ^ 0, when we can hardly find the analytical dependence of 
t upon T], the use of the formula ( |2.3| ) can help to find the dependence of t upon r] by 
quadratures: 



t = const + 



dp 
P 



^(u/P) + {2B/(3cl)p^ + {cl/cl)p\ p=ic, + C2v) 



Wot a = 1 it seems more convenient to substitute the representation ([4.3|) directly into Eq. (p.2|) 
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and formulas ( p^.ll ), ( 2.12 ) have the form 



« = {c2/p)^J{u/(3) + {2B/f3cl)p^ + {cl/cl)p\ v = \n{p/a). 

The behavior of a and v upon t, resulting from Eqs. ( [4 .61 ) for i? = as well as 
(|2.11| ), (|2.12|) for i? 7^ will be obtained by numerical integration. Some characteristic 
plots, related to these equations and ( [4.7| ), are shown in Fig. ||. 

(b) In the general case when i? 7^ we use the following representation for g, p and s 
00 00 ^00 

g = J]A„M^ p = ^S„M^ s = so--ln(^C„M"), (4.8) 

n=0 n=0 ?i=0 



while the remaining dependencies are obtained by integrating the equations in ( |4.1| ) 
for t and v. Substitution of (^^ ) into ( ^!1D yields an infinite set of equations for the 
coefficients: 



j=0 



2(1 + u)6k,3 + f35k,i - PBk - Be-^^'C, 



k 



1 — a 



^ ((fc + 1 - j)Bu+i-,Aj + 25fc_,5,) 



(4.9) 



j=o 



^ (4,2 + A + 5fc_i(l-5fc,o)) =0, 



^ ((fc + 1 - ])Ck+i-jAj + 3Cfc_,5,) = . 



Again, we present two examples of solutions that correspond to truncated series ( |4.8|) 
for some k > fcmax- 

(bl) The first solution corresponds to a = 0, A; max = 2 

Ao = Bo = Co = Ai = Ci=B2 = C2 = 0, A^ = Bk = Ck+i = , k>3, 
5i = l, A2 = -l, C3 = l, 

s = -{l/3)\n{-7ul/9B) + \n{l + uoit-to)) , v = Vq = const. , (4.10) 

Mo Uo 
u = p = r , q 



l+Mo(t-to)' \l+Uo{t-tQ) 

Here again we see the effect of explosive behavior for mq < at t — to ~(l/'^o) a^nd 
asymptotic stability at t — *■ cxo in the opposite limit Mo > similar to the observation 
above. Again, the curvature is positive and approaches 
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Figure 5: Characteristic behavior of a{t) and v{t) resulting for B = from (a) equations 
(^) for C = 1 and Vo = I (first, top panel), (b) equations Q for = 100, G = 1 
(second panel); and from equations (|2.11|J2.1^ ) for = 100, C2 = 1, Cg = 0.1 with (c) 
i? = 10 (third panel) and (d) B = —0.45, (bottom, fourth panel), respectively. 
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(b2) The second solution corresponds to some a 7^ 0, fcmax = 2: 

Ao = 5o = Co = = 52 = C2 = , A, = 5,, = Cfc = , k>3, 



Ci = ^(l-i?i)exp(3so) , «=(^— ^J=l + 4 , 
s = -(1/3) ln(/5(l - Si) + ^ In (1 + 'iB^u^it - to)) , (^-H) 

v = VqA — — — In (1 + 3i?ino(t - to)) , ^ 



35i ' ' ^" ' 1 + 35iUo(t - to) ' 

^ _oo / !fo Y % 

^ '^l + 3i?iizo(t-to); ' ^ H + 35i«o(t - to) ■ 

Here, the system demonstrates explosive behavior for B^uq < 0, while for BiUq > 
we have vanishing u, q and p for t — 00, and the raise of s and f is softer than in the 
stationary case. Here, the curvature is negative and approaches — (2/3)t~^. 

4.2 Series in q 

Here we will examine only the particular case of 5 = 0. Then we use the following 
representation for u and p 



00 00 

u = 

n=0 n=0 



X^^ng", p = J2Bnq\ (4.12) 



while the remaining dependencies are obtained by integrating the first three equations 
in ( |4.2| ) for t, s and v. Substitution of the expansions (|4.12|) into the next two equations 



in (^4.2|) yields the following infinite set of equations for the coefficients: 

k 

+ 1 - 3)Ak+i-jRj - 6k,i = , 

i=o 

k 

+ 3)Bk+i-jRj + '^Bk-jB^ (4.13) 



3=0 

1 / A 

0, 



1 ^ 

j=0 
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where Rj = (5,, o - 1) (^35,_i + j + - Aj) 

i=0 i=0 ^ ^ 



(4.14) 



Despite being of principal interest this set of equations is quite comphcated. There 
exists one evident truncation of the series (|4.12| ) namely when only zero order terms 



are taken into account, Aq 7^ 0, i?o 7^ 0, whilst A^ = = for k > 1. This truncation 
corresponds to g = and the stationary solution (|3.8| ). Unfortunately, in the case g 7^ 
we were not able to find a nontrivial truncation of these series. 



5 Stability analysis 

5.1 Small perturbations 

Here, we consider only the case B = 0. Then the basic equations have stable solutions 
of the form (|3.6|) . By linearizing these equations with respect to small perturbations 
6u = u — u St, = Q, = P — P St we obtain a system of linear equations that admit 
solutions of the form 6u (x U e^*, 5g oc Q e'^* and 6p oc P e'^* where A obeys the 
characteristic equation 

A' + A f 3p ,t + + T^{2p St + 3ust)+P- 6(1 + u)u\, 

\ 1 — a / I — a 

X f A + 4p 3t - ^] - ^ (p - ^] (A + p 3t + 2w St) = . (5.1) 
\ 1 — aj 1 — a \ I — a J 

Below, on Fig. H, we present a plot illustrating the dependence of real parts of solutions 
of the equation ( p.l| ) upon the parameter uj. The graphic shows that there exist stable 
solutions in a large region of the parameter u. 

The corresponding numerical solutions of input equations (p.5|) on Fig. |^ shows the 
behavior of the functions s and v in the vicinity of the stationary state. Straight lines 
(curves 1) correspond to stationary solutions (|3.6| ) with m st = 3, 19549, p st = 3, 20687 
for the left panel and Ust = —138, 843, p st = 69, 5447 for the right panel. For the left 
panel curves 2 relate to — st = —0.02, po — P st = 0.01 and go = 0.01; curves 3 
relate to Uq — u^t = —0.02, Pq — p st = —0.01 and go = —0.01. For the right panel 
curves 2 relate to uq — u = —0.02, po ^ Pst = —0.01 and go = 0.01; curves 3 relate to 
Mo - Ust = 0.01, Po - P st = -0.01 and go = -0.01. 

5.2 Transition to the stationary solution 

It is evident that for initial values that are close enough to the stable stationary values 
we can get an analytical solution for arbitrary initial values u{t = 0) = Uq, p{t = 0) = po 
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Figure 6: Dependence of real parts of solutions of the characteristic equation ( |0| ) upon 
u) that demonstrates the existence of stable (curves 1, 2 and 3 on the left panel corre- 
spond to negative real parts) and unstable (the curve 3 on the right panel corresponds 
to the positive real part) stationary solutions. 

and q{t = 0) = go- This solution is given by linear combinations of exponential 
functions defined by solutions of the characteristic equations presented in the previous 
subsection with coefficients, depending upon the initial values 

q = Di exp(Ait) + D2 exp(A2t) + -D3 exp(A3t) , 

u = u,, + (/^i/Ai) exp(Ait) + {D2/\2) exp(A2t) + [D^/X^) exp(A3t) , 

p = Pst + DiAi exp(Ait) + D2A2 exp(A2t) + D-sAs expiXst) , (5.2) 

. {pst + 2ust + \i) / Xi 

with Ai = — . 

(4pst + Xi){l - a) - 

Here the constants Di, D2 and D3 are expressed in terms of the initial values as follows: 

D. = ^, fl, = ^. = (5.3) 

ZAo ZAq ^^0 
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Figure 7: Dependence of functions s and f for u; = 500, = 10 for stable (left panel) 
and unstable (right panel) solutions. 



Integrating the remaining equations we get the following formulas for s and v that 
describe the transition to the stationary state 




It is clearly seen that the asymptotic behaviour of the functions s and v depends linear 
upon t, and the slope of these lines is defined by the stationary values, while the slope 
at the point t = is defined by the initial values. 

In order to compare the perturbative solutions ( |5.4| ) with the exact ones of the basic 
equations ( |3.5| ) we considered the behaviour of s, v, u, p, q and a for initial values uq, 
Po and go that are close to the stationary values Ust, Pst and g = for a; = 500, B = 



20 



and different values of A^. An example, for = 10, is presented in Fig. |^. 

These plots demonstrate that there is only a small discrepancy between exact and 
approximate solutions for physically interesting quantities a (resp. s) and f , whereas 
the difference for m, p and q is somewhat larger. These differences decrease for increasing 
values of A^. Obviously, the exact solution is approximated very fast. 

But, even if the initial values differ substantially from the stationary ones the 
approximate formulas still give a satisfactory accuracy for uq = Pq. Moreover, as the 
input equations are linear in q it appears that the approximate solution is sometimes 
valid even when the initial value of go is not small. This is clearly seen in Fig. |^. 

Also in these cases the stationary solution is approximated after a few units of time. 

6 Conclusion 

Analysing the evolution equations in conformal time we were able to present a new 
exact solution (for a = 1 ~ = 7/6). After transforming to physical time we found 
'stationary' solutions of the corresponding basic set of first order differential equations 
(restricted to -B = 0) which generalize the inflationary solution already known from an 
earlier study . The numerical analysis of the basic set showed quite different types of 
dilaton-driven evolution in cases B = as well as -B 7^ 0: In a wide range of the param- 
eters of the theory the stationary solution is approximated very fast; furthermore, there 
are also quite different non-stationary solutions showing explosion type, oscillatory as 
well as collapsing behaviour. Seeking solutions as (finite) power series we were able 
to find additional particular solutions of the basic set with a (modified) logarithmic 
raise of s{t) which, in comparison with the stationary solution raising linearly, lead to 
a more soft inflation; however, nontrivial dilatonic behaviour results only in the cases 
(a2) and (b2). Finally, we presented a stability analysis for the case B = thereby 
showing how the transition to the stationary solution occurs. 

It is very interesting that the mathematical methods applied here to study the 
conformally fiat Universe with time-dependent dilaton solutions leading to effective 
forth order equations of motion are general enough to be used also in various related 
contexts. As an immediate application we may consider the problem of quantum 
creation and annihilation of Anti de Sitter Universe and Anti de Sitter black holes 
due to quantum effects of (dilaton coupled) matter. Such a phenomenon has been 
investigated recently in Ref. |^. In these works the anomaly induced effective action 
for dilaton coupled spinors has been used similarly to the present paper. Moreover, 
the Anti de Sitter Universe in conformally fiat coordinates looks very similar to the de 
Sitter (inflationary) Universe: the corresponding equations of motion are - with some 
change of signs - almost the same as Eqs. ( ^.6] ) and ( p^.T] ) (where, however, the role of 
time is played by the radial coordinate in AdS Universe). Hence, our results may be 
used in the construction of other variants of the asymptotically AdS solutions which 
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Figure 8: Exact numerical (solid line) and approximate analytical (|5.4]) (dashed lines) 
solutions of basic equations ( p.5|) for oj = 500, N = 10, B = 0, qq = 1 and initial values 
■"o = Po = 2 close to stationary ones Ugt = 3.19549, Pst = 3.20687. 
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Figure 9: Exact numerical (solid line) and approximate analytical (|5.4]) (dashed lines) 
solutions of basic equations (|3l5| ) for = 500, = 10, -B = 0, go = 50 and initial 
values uq = po = 0.5 not close to stationary ones Ust = 3.19549, pst = 3.20687. 
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have been found in Ref. 0. 
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